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Abstract 
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scribed by the 7 metric. Their behaviour is compared with the spher- 
ically symmetric situation, bringing out the sensitivity of the trajec- 
tories to deviations from spherical symmetry. 
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1 Introduction 



The influence of small perturbations of a Schwarzschild black hole, on the 
trajectories of test particles around the source, has atracted the attention 
of many researchers. These perturbations are usually introduced as addi- 
tional mass and charge concentrations |Ij, || 0, |], j|, H, magnetic fields [[F| or 
gravitational waves ||. The common results of all these works being, that 
essentially any perturbation of a Schwarzschild black hole will lead to chaotic 
orbits. 

In this work we want to present another approach to the problem of intro- 
ducing perturbations to spherical symmetry. This consists in considering an 
exact solution of Einstein equations continuously linked to the Schwarzschild 
metric, through one of its parameters. The rationale behind this approach 
lies on the known ||, though usually overlooked, fact that as the source ap- 
proaches the horizon, any finite perturbation of the Schwarzschild spacetime, 
becomes fundamentally different from any Weyl metric, even if the latter is 
characterized by parameters whose values are arbitrarily close to those cor- 
responding spherical symmetry. 

The solution to be considered here is the so called 7 metric [|H], |TTJ, which 
is also known as the Zipoy-Vorhees metric [12 , and belongs to the family of 
the Weyl solutions p| . 

The motivation for this choice may be found in the fact that the 7 metric 
corresponds to a solution of the Laplace equation, in cylindrical coordinates, 
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with the same Newtonian source image [II , as the Schwarzschild metric (a 



rod). In this sense the 7 metric appears as the natural generalization of the 
Schwarzschild spacetime to the static axisymmetric case. 

We shall find the geodesic equations for test particles in the 7 metric. 
Particular attention will be devoted to circular and radial geodesies. The 
qualitative differences in the dynamics of the test particles as compared to 
the spherically symmetric case will be illustrated and discussed. 

The paper is organized as follows. In the next section the 7 metric is 
briefly presented. In section 3 geodesic equations are found and analyzed 
and in section 4 the gravitational and centrifugal forces are studied. Finally 
the results are discussed in the last section. 

2 The 7 metric 

In cylindrical coordinates, static axisymmetric solutions to the Einstein equa- 



tions are given by the Weyl metric [I3~| 



with 



and 



ds 1 = e 2X dt 2 - e" 2 VW + dz 2 ) + p 2 d<p 2 }, (1) 



A,pp + p' l \ p + X, zz = 0, (2) 



V, P = p{^ 2 p ~ ^ 2 J> = 2pA iP A jZ . (3) 

Observe that (2) is just the Laplace equation for A in the Euclidean space. 



The 7 metric is defined by [|Kj 



-,2A 



,2/i 



R 1 + R 2 -2my 

(gi + R 2 + 2m)(R 1 + R 2 - 2m) 
4,R±R 2 



where 



R\ = p 2 + [z - m) 



R\ = p 2 + (z + m) 



(4) 
(5) 

(6) 



It is worth noticing that A, as given by (4), corresponds to the Newtonian 
potential of a line segment of mass density 7/2 and length 2m, symmetrically 
distributed along the z axis. The particular case 7 = 1, corresponds to the 
Schwarzschild metric. 



It will be useful to work in Erez- Rosen coordinates |L2|, given by 



p 2 = (r 2 — 2mr) sin 2 9, z = (r — m) cos#, 



(7) 



which yields the line element (1) with (4,5) as [HI 



ds 2 = Fdt 2 - F'^Gdr 2 + Hd6 2 + (r 2 - 2mr) sin 



2 /ijj.2 



(8) 



where 



F 
G 
H 



;| 2m\ 7 



r / 
^.2 



2mr 



r»-i 



2mr + m 2 sin # 



2mr)' 



(r 2 — 2mr + m 2 sin ^) 7 1 



(9) 
(10) 

(11) 
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Now it is easy to check that 7 = 1 corresponds to the Schwarzschild metric. 
The total mass M of the source is HID, O 



M = 7m, (12) 

and the quadrupole moment Q is given by 

g = 2(l_ 7 2 )M 3. (13) 

From (13) we have that 7 > 1(7 < 1) corresponds to the oblate (prolate) 
spheroid. 

3 The geodesies 

The equations governing the geodesies can be derived from the Lagrangian 

2£ = g aP x a x^ (14) 

where the dot denotes differentiation with respect to an affine parameter s , 
which for timelike geodesies coincides with the proper time. Then, from the 
Euler-Lagrange equations, 

d_ (dC \ _ 9£ 

ds \dx a J dx a 

it follows, for the metric (8), 



0, (15) 



i + jtri/ = 0, (16) 
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-2r 



. 9 



A? 1 



7 



e-vi/ 



At 2 



r + (l-7 2 )m 2 ^r^(- + ^)sin 2 ^ 



-4r0(l - -f 2 )m 2 



e -(l+7)" 



+0 2 



r 2 e (l-7)^ f2 



+ (1- 7 K 



r 2^ 

+ (l- 7 2 )m 2 
2 



sinflcosfl- tV'e 7iy 
+ v sm 2 6 1 



A7 2 



+ ^2 r 2 e (l- 7 ). 



+ (1-1)1/ 



sin 2 = 0, (17) 



+6»f 



_ „2 (1-7)1/ _ -71/ 

-e —-^- - 2 (1 - 7 2 )m 2 sin cos 6 



r 2 e (l-7)" 



A 7 2 -l 

+r 2 (l — 7 2 )m 



- + (i- 7 y 

Lr 



+ (1 - 7 2 )m 2 ^ + v< ) sin 2 



e -(l+7)f 

- ■ :: sinflcos# + 2 r 2 e {1 - 7)iy sin#cosfl = 0, (18) 



with 



+ (1-7)1/ +2^cot^|e (1 - 7)iy sin 2 = 0, (19) 



2m 



A= 1 + 



m 2 e v sin 2 6> 



(20) 



and prime denotes differentiation with respect to r.It is a simple matter to 
check that if 7 = 1 we recover the geodesic equations of the Schwarzschild 
spacetime. 

Let us first consider circular geodesies. From (16-19) we obtain using 
r = 6 = 0, 



t = <f> = 0, 
7 tW = <j) 2 r 2 e (1 ^> 



- + (i- 7 y 

r 



(21) 

sin 2 6, (22) 
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2<j) 2 r 2 e (1 -~ f)u sm9cos9 = 0. (23) 

Then from the defmiton of the angular velocity w of a test particle along 
circular geodesies 



" = f , (24) 



we obtain, using (22), 



2 7m / 2m\ 27 1 



^ = 1 - — . (25) 

[r - (1 + 7)m]r 2 sin 9 \ r J 

First of all, we observe that, as it follows from (23), circular geodesies out of 

the equatorial plane 9 = tt/2, are now possible, if only 7 < 1 and r = 2m. 

This kind of trajectories do not exist in Schwarzschild spacetime (7 = 1). 

From (25) it also follows that physically meaningful values of uj at r = 2m, 

exist only for 7 = 1/2. Therefore circular geodesies for 9 7^ 7r/2 can ocur if 

only 7 = 1/2 and r = 2m, with angular velocity 

^ 2 = 9 „ . (26) 

4m 2 sin 2 9 

Circular geodesies outside the equatorial plane also exist in the Kerr metric 
15]. These kind of orbits implies the presence of repulsive forces whose 



nature is still not well understood 16, 171. However it should be stressed 



that since r = 2m represents a physical singularity in the 7 metric, these 
orbits are deprived of physical meaning. In the weak field limit one obtains 
from (25), for 9 = tt/2, up to the first order in m/r, 



r, 7m 



l + (l- 7 ) — 

r 



7 



3 "' Of™)*. (27) 



r 



If 7 = 1, we recover the well known Kepler law, which is also valid, if 7 = 1, 
without any approximation. 

If 7 7^ 1 the second term within the brackets, in (27), gives the correction 
due to the quadrupole moments of the source. In troducing the dimensionless 
variable x = 2m/r, we can rewrite (25), for 9 = 7r/2, 

uj 2 = m 2 uo 2 = A[/ n . — — ; -. 28 

4[(l- 7 )x + 2(l-x)] V ; 

Figure [l] shows the behaviour of uj 2 as a function of x for different values of 

7, including 7 = 1. 

The bifurcation between the 7 = 1 and 7 ^ 1 cases, for values of x close 
to one, clearly illustrates the sensitivity of the system under perturbations 
of 7 in the neighborhood of 7 = 1. 

For comparative purposes it will be useful to find an expression for the 
tangential velocity of the test particle along the circular geodesic. From 
flUI , we have 

W a ^ 29 
y/goodx ' 

with 

V a = (0,dx\dx 2 ,dx 3 ), (30) 
and we obtain from (8) for 9 = n/2, 

W 2 ee W a W° = (l - ^ u 2 r 2 = (31) 

V t J r — (1 + 7)777. 

where uo 2 is given by (25). In the weak field limit, m/r 1, the classical 

expression, W = ur is recovered. Figure ^ shows W 2 as a function of x 



for different values of 7. Furthermore, we see from (31), that null circular 
geodesies appear when r = (1 + 27)771. Since the physical singularity lies for 
r = 2777 then null circular geodesies can only exist for 7 > 1/2, otherwise, for 
7 < 1/2, they do not exist. 

The Levi-Civita metric, which represents the field of an infinite line mass 
of constant energy density a, can be obtained as a limiting case from the 7 
metric by taking m — > 00 and associating 7 = 2a |TjJ. Next we see, in the 
equatorial plane, how the angular velocity u and the tangential velocity W 
of the circular geodesies, respectively given by (25) and (31), behave for big 
values of m. 

Rewriting (25) in terms of the coordinates given by (7), with 6 = ir/2, 
we obtain 

JL = 7777(y/777 2 + fj 2 - m) 2 ^ 1 

{^/mF^T]? — 7777) (777 + y/rri 2 + p 2 ) 2 T +1 
Expanding in series (32) for big values of m and keeping only the term in its 
lowest order of p/777, we obtain 

2a 



UJ 



2 - -^—(AmY^P 2 ^^, (33) 



2a 



where we have substituted 7 = 2a. The expression (33) is the same as for 
the Levi-Civita metric obtained in [[[]]] provided the topological defect a, as 
given in [ljj, is associated in (33) as a = (4m)~ 2cr . In [[0| it is proved that 



the limiting case of Levi-Civita spacetime from the 7 spacetime produces an 
infinite topological defect. 



9 



Now rewriting (31) in terms of the coordinates given by (7), we have 

W = . (34) 

\Jmr + fr — 7m 

Expanding (34) in series for big values of m up to the second order of p/m, 
we obtain 



where we have substituted 7 = 2a. From (35) we have that, for a given 7 
and a fixed radius p, decreasing the length 2m of the line segment of mass 
decreases the tangential speed W of the circular geodesies. When m — > 00 we 
have that the tangential velocity (34) is the same as the one obtained for the 
Levi-Civita metric [ 17] , being 



lim w 2 = = (36) 

m->oo 1 — 7 1 — 2(T 

Observe that (36) sets a constraint on possible values of o to avoid values 
of W larger than 1, i.e. the velocity of light. This is in contrast with the 
situation in the 7 metric (34), where such constraint W < 1 involves 7 and 
the radial coordinate p of the orbit. 

Let us now consider radial geodesies. From (8), in the 9 = it/2 plane, it 
follows 

. -.2 1 

9 / 2m \ 7 9 / 2m \ 7 / r — 2rmr 
2 ' 1 - r 2 1 



1 = r 1 



r / V r / \r 2 — 2mr + m 2 



_^ (1 _^y-v. ,37) 



r 
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Next, it follows from (15) 
dC 



2m\ 7 



constant — E = t ( 1 ) . 

dt V r J 

dC ■ / 2m\^ 7 2 

— - = constant = L = —<p 1 r , 

d<f) \ r J 



(38) 
(39) 



where E and L represent, respectively, the total energy and the angular 
momentum of the test particle. Then using (38,39) in (37), we obtain 

1 (E 



- {E*-V% 



(40) 



where V(x), which can be associated to the potential energy of the test 
particle, is given by 



V 2 



_ / 2m\ 7 
r / 



m 2 L 2 ( 2m\ 7 ~ 1 



(41) 



r 2 — — 

— 2' 



(42) 



with 

or, in terms of x (41) becomes 

v = (i _ (i _ ^r 1 + ii . ,43, 

V xj \_x A \ xj 

Figure |3] shows V for different values of 7 and L 2 = 55. For indicated values 
of the total energy E, there exist unstable circular orbits, which illustrate 
the sensitivity to small perturbations. 
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4 Gravitational and centrifugal forces 

The study of the gravitational and centrifugal forces follows the same treat- 
ment as in |[20|| . The four velocity of a particle moving on a circular orbit in 
the equatorial plane 6 = n/2 of the 7 spacetime (8) can be expressed as 

«« = (tf + fi<5J)I\ (44) 



where Q is a parameter and, since v a v a = 1, T is given by 



(45) 



F 2 - fi 2 r(r - 2m)' 

Assuming uniform motion, Q =constant, the acceleration a a of the particle, 

a a = iA,« , (46) 

using (8,44,45) is 



x 



r -2 7 ( r _ 2m) 2 f- 1 - Q 2 r 
2m\ 2 ^ 1 



T1U-!!!1\ - n »[ r -(l +7H . (47) 



7™ ^ V T* J 



If 7 = 1, then (47) reduces to the Schwarzschild spacetime result pO 



^1 / 2m \ m — fi 2 r 3 . , 

a = S r~ 1 5 ™-s ■ 48 

r V r J r — 2m — irr 6 

The acceleration a" can be split into its gravitational component g a , when 
f2 = 0, and its centrifugal component 

-a a = 9 a + c a , (49) 
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and so we have from (47), 

9 = - 5 r-^{ 1 -—) {}--) > ( 50 ) 

When the gravitational acceleration is balanced by the centrifugal accelera- 
tion, which means a a = 0, then we have Q = uj, as it can be checked from 
(25). The gravitational acceleration a a always points inwards, being g a < 0. 
The centrifugal acceleration c a for values of r > (1 + 27)771 is always positive, 
while for r < (1 + 27)777, it becomes negative. Since r = 2m is the physical 
singularity only 1/2 < 7 allows negative values for c a , while 1/2 > 7 makes 
always c a positive. 

5 Conclusions 

We have considered deviations from spherical symmetry by considering exact 
solutions of the Weyl family, instead of perturbing the Schwarzschild metric. 



As it should be expected from the Israel theorem pl| , these approaches 
differ qualitatively as the orbit of the test particle gets close to the horizon. 
Though it was not our purpose here to describe the chaotic behaviour of 
the trajectories, it should be clear that such behaviour is the expression 
of the sensitivity to small changes of 7. This sensitivity in turn, appears 
sistematically in the cinematics of the particles, for orbits close to 2m, as 
illustrated by figures |l]-|3|. 
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Figure 3: V as function of x for three different values of 7. For the three 
values of E, there are three different unstable circular orbits. 
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